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Abstract. The statistics of the average height fluctuation of the one-dimensional 
Kardar-Parisi-Zhang(KPZ)-type surface is investigated. Guided by the idea of local 
stationarity, we derive the scaling form of the characteristic function in the early-time 
regime, t <C TV 3 / 2 with t time and N the system size, from the known characteristic 
function in the stationary state (t ^> TV 3 / 2 ) of the single-step model derivable from a 
Bethe Ansatz solution, and thereby find the scaling properties of the cumulants and 
the large deviation function in the early-time regime. These results, combined with 
the scaling analysis of the KPZ equation, imply the existence of the universal scaling 
functions for the cumulants and an universal large deviation function. The analytic 
predictions are supported by the simulation results for two different models. 

Keywords: Kinetic growth processes (Theory), Fluctuations (theory), Integrable 
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1. Introduction 

The Kardar-Parisi-Zhang (KPZ) equation pj is a nonlinear differential equation with 
noise representing in the simplest way the roughening phenomena of growing surfaces. 
The roughness is measured by the spatial fluctuation of the surface profile, and an 
initially flat surface of finite size finds its width increase in the early-time regime and 
then saturate in the stationary state. The KPZ equation offers a theoretical framework 
to understand such a kinetic roughening, and many real-world nonequilibrium interfaces 
as well as computational models belong to the KPZ class or its variants [2]. 

A remarkable universal quantity for the (l+l)-dimensional growing surfaces is the 
large deviation function (LDF) [SJEJElin]- It concerns the tail behavior of the probability 
distribution of the spatially-averaged height, h = iV _1 Y*£=i hi with iV the system size 
and hi the height at a position i. In the stationary state, it has been claimed to be 
universal within the KPZ class. The exact LDF [3] derived for the single-step (SS) 
model [Zj has the cumulants' ratio (/i 3 )c/((^ 2 }c(^ 4 )c) = 0.41517..., which has been 
reproduced for other computational models with no significant deviations jU E|. The 
fluctuation of the average height from configuration to configuration does not saturate 
but all of its cumulants grow with time even in the stationary state and the LDF has 
asymmetric tails. The LDF of the average height in the stationary state is of interest 
also in the theory of nonequilibrium systems. It has been shown that the probability 
distribution of the phase space contraction (or entropy production) rate has a non- 
trivial symmetry property in dissipative chaotic systems and this has been considered 
as a generalized fluctuation theorem for nonequilibrium systems 0. Recent studies 
show that the theorem is applicable also for stochastic dynamics of interacting particles 
with the symmetry property found in the LDF of the particle current, which corresponds 
to the total height increase, Nh, in the language of the growing surfaces. 

The LDF of an averaged quantity contains a lot of information specific for 
a given system. For example, one can obtain the dynamic correlations among 
individual variables, which are not revealed by the Gaussian behavior around the mean 
value. Previous works, however, have been restricted to stationary states, where, in 
mathematical terms, only the ground state of the time-evolution operator survives and 
thus its LDF could be derived analytically. In this work, we explore the LDF of the 
average height in the early-time regime of the KPZ surface. By the early-time regime, we 
mean N — > oo, t — > oo, and t <C N z , with z = 3/2 the dynamic exponent. In the early- 
time regime, a flat surface evolves into its stationary state which is rough. The LDF of 
the average height corresponds to the Legendre transformation of the logarithm of the 
characteristic function, as the entropy to the free energy. We derive the characteristic 
function in the early-time regime of the SS model from its stationary-state counterpart 
that is exactly known [31 EJ EJ - To do so, we use the assumption that the scaling 
property and the model-parameter dependence of the latter are the same as those of 
the locally-stationary segments in the early-time regime, the size of each of which scales 
with time as t x l z . Combined with this result, the scaling analysis of the KPZ equation 
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enables us to see how the characteristic function depends on the parameters of the KPZ 
equation. We derive the LDF and the dynamic scaling behavior of the cumulants from 
this characteristic function, and their universal parts are identified in the numerical 
results from two different models in the KPZ class. 

2. LDF in the stationary state of the SS model 

The (l+l)-dimensional KPZ equation describes the evolution of the heights h(x,t) on 
a substrate {x\0 < x < L}: 



where the noise r}(x,t) satisfies (rf(x,t)r)(x , ,tf)) = 2D5(x - x')5(t - f) p. It is well 
known that the rescaling of the space x — > bx gives rise to h — > b a h and t — > b z t with 
a = 1/2 and z = 3/2 as a result of the Galilean invariance and the fluctuation-dissipation 
theorem valid in one space dimension [TUj . 

We first review known results about the LDF derived from the Bethe Ansatz 
solution. The SS model [Zj is a discrete growth model belonging to the KPZ class, 
where the height hi for 1 < i < N can take even (odd) integer numbers for i even (odd) 
and satisfy the SS condition, \h i+ \ — hi\ = 1, for all i with /i/v+i = hi. The height at a 
randomly chosen site is increased by 2 with probability (1 + e)/2 or decreased by 2 with 
probability (1 — e)/2, with < e < 1, provided the new height does not violate the SS 
condition, and the time t increases by 1 after N such attempts. This model is equivalent 
to the partially asymmetric exclusion process. Hence the time-evolution operator for 
the height step configuration {cr} with o~i = hi — hi-i is identical to the asymmetric XXZ 
spin chain Hamiltonian which in turn is diagonalized by the Bethe Ansatz [TT) IT2"] . The 
ground state energy of a parameterized version of the Hamiltonian then determines the 
characteristic function of the total height increase, J = J2i(hi(t)—hi(0)), for tN~ 3 / 2 ^> 1. 
Note that the excitation gap scales as iV -3 / 2 (z = 3/2) [IUH2]- With (. . .) denoting the 
expectation with respect to the distribution of J, P(J), the characteristic function of J 
is given by [3J HJ IS] 



for 7iV 3 / 2 finite, where (J) = 2ep(l — p)N 2 t(N — 1) _1 and the function G(x) is given 
parametrically as G = Y^=i{—z) n {n^l 2 —n~ b l 2 ) andx = — Y^=\{— z) n n~ 3 / 2 for x > — 1 
along with its analytic continuation for the region x < — 1 r 5j. From now on, hi(t) refers 
to hi(t) — hi(0) and h = J/N. Then the characteristic function of the average height 
distribution is simply, (e^ h ) = e <h)+g(y/N,N,t) with ^ = (j)/jv. 




(1) 




el{J)+ g (l ,N,t) with 
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The logarithm of the average height distribution, In P(h), is related to ln(e 7/l ) via 
the Legendre transformation: 

In P(h) = ln(e 7/t ) -7/1 and h = ^- ln(e 7 ^). (3) 

Using equation (0), one finds that for t ^> N 3 ^ 2 , 
P(h = (h) + yN/t) ~ e~ ts(y) with 



J_ ^l-„) / y \ 
ny > JV'/ 2 2^ \2cp(\-p))' w 

where the function H(x) is defined parametrically as H = [G(z) — z(d/dz)G(z)]/\/2 
and x = (d/dz)G(z) and behaves as H(x) ~ |x| 3 / 2 for x — > —00 and H(x) ~ x 5//2 for 
x — > 00 j^j. 

3. LDF in the early-time regime of the SS model 

3.1. Local stationarity 

A feature of the stationary state of the one-dimensional KPZ surface is that the height- 
steps at different sites and at equal or different times are uncorrelated as the fluctuation- 
dissipation theorem states [TX3 - On the other hand, the height increases, h(t + At) — h(t) 
with different t's may be correlated. For example, a flat surface has many sites where 
the heights can increase while in a V^-shaped surface the heights has only a few such 
sites. Moreover, such morphological features relevant to the height increase can survive 
over a certain period of time £t. The factor t/N 3 ^ 2 in Q{p/,N,t) implies that ~ N 3 ^ 2 
[figure l^a)]. This motivates us to decompose the total height increase J into t/N 3 ' 2 
effectively-independent partial height increases, as J = Ya=o -^jv 3 / 2 (N 3 ^ 2 , N), where 

J t>x {At,Ax) = Y?[hi(t + Ai)-hi(t)]. (5) 

i=x+l 

In this notation, J = J Q0 (t,N). That such partial height increases are statistically 
independent is our main assumption. To check its validity we calculate, for the SS 
model, the normalized correlation defined as 

= (^(At,*QAWAt,*))c (6) 
^(J tfi (At,N) 2 ) c (J t+Atfi (At,N) 2 ) c 

around At ~ N 3 ' 2 where (AB) C = ((A - (A))(B - (B))). The results shown in 
figure IHb) exhibit a rapid decay of the correlation for At ^> N 3 ^ 2 and confirm the 
mutual independence of J iN 3/2^(N 3 ^ 2 , N) with different z's. 

In the early-time regime when the roughening is still ongoing, the height-steps at 
time t are uncorrelated only within each segment of size £ x ~ t 2//3 , and such locally- 
stationary segments expand as time goes on and merge into one in the stationary state. 
At a given time t\ with 1 ti iV 3 / 2 , we may assume that most of the height 
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Figure 1. Correlation of the partial height increases, (a) Example of a stationary 
correlated domain in the stationary state, D, and that in the early-time regime, D, 
on both of which all the partial height increases are correlated and the height-steps 
are uncorrelated. (b) Plot of Cti me (Ai) versus AtN~ 3 / 2 for the SS model with e = 1 
and p — 1/2, and the system size N = 20,40, and 80. (c) Plot of C S p ace (Aa;) versus 
Axt- 2 / 3 for the SS model with e = l,p= 1/2, N = 10 4 at time t = 50, 100, and 200. 



steps are uncorrelated and all of the partial height increases are correlated within a 
domain [0, t\] x [x,x + £ x ] with £ x ~ t^ 3 (D in figure ^a)), in the same way as within a 
stationary correlated domain [t 2 , £2 + £t] x [0, N] of the stationary state with £ t ~ N 3 ^ 2 
and t 2 ^> N 3 / 2 (D in figure ^a)). To check this idea, we measure for the SS model the 
normalized correlation function, 

,x+Ax (t,Ax)) c _ 

^ space \ ) — j 1 \ < ) 

(J 0>x {t, Ax) 2 ) c (J 0iX+ Ax(t, Ax) 2 ) c 
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as shown in figure dc). The decaying behavior of C space (Aa;) confirms the mutual 
independence of the partial height increases, Jo^it, £ x ) an d Jo,j^ x (t,^x) on different 
(i 7^ j) domains of size tx^. 

3.2. Dynamic scaling and LDF in the early-time regime of the SS model 

The same correlation landscape between the stationary correlated domains of the 
stationary state and of the early-time regime implies that the characteristic function 
in the stationary state, e 7 ^ + ^ 7,Ar '^ gives the characteristic function for a stationary 
correlated domain of the early-time regime, ^e" ,Jo - x ^ t ^ x ' ) ) , provided the system size N 
in the former is replaced by the correlation length £ x ~ at 2 / 3 in the latter with a a 
constant. In this mapping, actual exact functional form of G(x) in equation (2) may 
possibly change. Then the characteristic function for the early-time height increase is 
determined as 

(e 7j ) = (e 7 ^^ J °>^ («*)) ~ ( e 7Jb,.(t,&))iv/& ^ 

We then have, with t < N 3/2 , 

(e 7j ) = e ^J)+G(y,N,t) with 



S(7^,t) = - x G (V27 7v /4p(l -p)~a^t) (9) 



for jt finite with a new function G(x) introduced. 

The expressions of Q{pf, N, t) and Q{f, N, t) together determine the dynamic scaling 
of the average height by the relation (h n ) c = N~ n (d n /d , y n ) ln(e 7j ) | 7= o, which yields for 
n — 2, 3, . . ., 

^ = tN in-m e[Ap{l _ p)](«^)/^„ {j^J^j With 

Here G( n )(0) = (5 n /^ n )G(^) \ z=0 and the same for & n \0). The scaled data of the 2nd, 
3rd, and 4th cumulants in the SS model with e = 1 and p = 1/2 are plotted in figure 121 
which agree with equation (fTUj) . 

From equation (jOJ) one also finds that the LDF for N ^> t 2//3 is given by 

P(h = (h) + y/t 1/3 ) ~ e~ Nfm with 
where is defined by H = [G(z) - z(d/dz)G(z)]/V2 and x = (d/dz)G(z). 
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Figure 2. Dynamic scaling of the (a) second, (b) third, and (c) fourth cumulants of 
the average height of the SS model with e = 1 and p = 1/2 (A = —1,A = 1) and 
the RSOS model [T3] (A = -0.75, A = 0.81) for N = 50 and 100. The data scaled 
according to equations (|TU)l and {T5)| are presented, showing good collapse and also 
confirming the small-x behavior of the scaling functions </>„ (x) . The legend is provided 
in (a). 
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4. Scaling analysis and universal LDF 

The LDF and the dynamic scaling behaviors shown in equations (fTT)|) and (jllj) can be 
represented in terms of the parameters of the KPZ equation which enables us to 
check their universality. Following the procedures proposed in Ref. [H], one finds the 
following dimensionless variables 

A, . \ 2 D \ 4 D 2 

h' = —h, x' = — -x, r = —t 12 

2u 2v 3 4z/ 5 v ; 

satisfy an equation without parameters, (d/dr)h! = (d 2 /dx 2 )h' + (dh/dx) 2 + rf{x',r) 

with (i]' (x[, Ti)r]' (x' 2 , r 2 )) = 8(x[ — x' 2 )5{ti — r 2 ). Therefore both Q and Q should be 

functions of jN(X/u) (for 7 J), (\ 2 D/u 3 )N, and (\ 4 D 2 /u 5 )t. Keeping the t and N 

dependence appearing in equations (j2J) and (JOJ), those functions turn out to be given by 

G(j,N,t) = |A|v / I-^i3(sign(A)v / I 7 iV 3/2 ) , (13) 



where A = D/u, and 



^^*) = p^7s^(AA7t). (14) 

Here, B(x) and B(x) are the scaling functions to be identified as follows. The nonlinear 
term with the parameter A in equation JIJ has been introduced to describe the lateral 
growth and thus A = (d 2 /du 2 )(d/dt)(h)\ u=0 with u the average slope of the surface [THj . 
which leads to A = — e for the SS model where (h) = 2ep(l — p)t and u = 1 — 2p. Also, 
the uncorrelated height-steps in the stationary state result in ((hi(t) — hj(t)) 2 ) c ~ A\i — 
j\ PHI, which gives A = 4p(l — p) from the comparison with the exact solution for the 
SS model [T7j. Substituting these values to equation (fT3*j) . one finds that equations (J2J) 
and ((TSJ) perfectly agree with each other under the relation B(x) = G(—\/2~kx)/ (2y / 27r). 
Also equations (0) and become equal under the relation a = \\\ 2 / 3 A 1 / 3 c with c an 
universal constant and the identity B(x) = G(— v / 27rxc 3 / 2 )/(2v / 27rc 5 / 2 ). 
With the functions Q and Q in equations (fT3*|) and (fUj). we find 

(^) c = tiV^^IAIA^ 1 )/ 2 ^ (lAIA 1 / 2 ^^ 2 ) . (15) 

To check the A and A dependence in equation (|TH|l and thereby explore the universal 
scaling function <p n (x), we performed numerical simulations of the RSOS model |13j . 
where the flat surface, hi = for all % at t = 0, grows via the integer- valued height at 
a randomly-chosen site being increased by 1 with probability 1 unless the new height 
differs from neighboring heights by larger than 1. The numerical estimates have shown 
that A ~ —0.75 and A ~ 0.81 for the RSOS model [TH]. Using these values, we plot the 
scaled data of the nth cumulants with n = 2,3, and 4 obtained from the RSOS model 
together with those from the SS model (e = 1 and p = 1/2) in figure El The good 
collapse of those data obtained for different times and system sizes strongly supports 
the existence of the universal scaling function <p n {x). 
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Figure 3. Universal LDF in the early-time regime. The data scaled according to 
equation (|17fl are shown, obtained from the simulations of the SS model (e = 1 and 
p = 1/2) and the RSOS model with N = 100 and 200 and t = 20 and 60. The 
collapsed data show the universal function (2y/nc 5 ^ 2 )^ 1 H (2cx) in equation fTTjl . To 
compare the functions H(x) and H{x), we draw the curve representing the function 
iJ(2xc)/(27r 1 / 2 c 5/2 ) with c = 1.2 chosen for the best fitting. 



The Legendre transformations of equations and (pPfj) give the LDF as 

for the stationary state and 

M = ^ | A |2/3 A i/3 2v ^g5/ 2 ^ (-2gsign(A) |A|1/ ^ 2/3 ) (17) 

for the early-time regime, respectively. While the acquisition of the data for the LDF of 
the stationary state is hampered by long evolution times, that of the early-time regime 
can be obtained within reasonable times. We plot \\\ 2 /' 3 A 1 / 3 t 2 /' 3 N~ 1 \n[P(h) / P((h))} 
versus | A | ~ 1/ ' 3 v4~ 2 / 3 t~ 1 / 3 (7z — (h)) with different system sizes and times in the SS 
model and the RSOS model in figure El which shows clearly a universal curve 
(2 % /ic 5 / 2 )- 1 ^(2£x). Note that A is negative for both models. 

Contrary to the stationary state, the exact functional form of the universal LDF 
of the early-time regime is not known yet since all the eigenstates of the time-evolution 
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operator are involved in the characteristic function. From the physics viewpoint, the 
detailed functional form of the probability distribution of the sum of correlated variables 
is expected to depend not only on the type of correlation but also on the type of the 
boundary condition 6J. Since the exact universal LDF of the stationary state has 
been obtained under the periodic boundary condition, it is not possible to draw the 
conclusion that H(x) = H(x) on the analytical basis. The boundary condition for the 
locally-stationary segment in the early-time regime is not obvious. On the one hand, 
however, the functions H(x) and H(x) are expected to display similar behaviors due 
to the same type of correlations among height increases and as seen in figure the 
function H(x) has indeed asymmetric tails as H(x) does. Moreover, we compare H(x) 
and H(x) in figure|3]by drawing the curve representing the function (2y / 7rc 5//2 ) _1 if(2a;c) 
with c = 1.2 chosen for the best fitting. The analytic curve and the numerical one are 
in good agreement over the whole range of data presented. This is consistent with a 
conjecture that H(x) = H(x). 

5. Conclusion 

We have investigated the statistics of the average height for the one-dimensional KPZ 
surface. The LDF in the stationary state implies that the average height is decomposed 
into statistically-independent components each of which is a set of correlated variables, 
and we checked its validity. To derive the LDF in the early-time regime, we assumed 
that the statistics of the average height in such a component remains the same across 
both time regimes, in that the scaling property and the model-parameter dependence of 
the characteristic function are not changed. The scaling analysis of the KPZ equation 
enables us to find the equation-parameter dependence of the LDF. Using this result 
in analyzing the simulation data obtained for two different models, we confirmed the 
universality of the LDF. 

We combined the exact results and the scaling analysis to explore the universality 
of the LDF, which can be complementary to the measurement of the cumulants' ratio. 
The exact functional form of the universal LDF is an open question. Given that little 
is known about the statistical properties of non- stationary systems, the present work 
could stimulate related research including the investigation of the fluctuations in non- 
stationary systems. 
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